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Dilepton emission rates from a hadronic gas at finite temperature and baryon density are com- 
pletely constrained by broken chiral symmetry in a density expansion. The rates can be expressed in 
terms of vacuum correlations which are measured in e'^e~ annihilation, r decays and photo-reactions 
on nucleons and nuclei Q,Qy{. In this paper, the theoretical results are summarized and the total 
dimuon yield is calculated by integrating the dimuon rates over the space-time history of a fireball 
based on hydrodynamic calculations with CERN SPS conditions. The resulting dimuon yield is in 
good agreement with the recent measurements reported by NA60 y|. 



I. INTRODUCTION 



o: 
o. 

P\| . It is expected that above a critical temperature, T^, w 170 MeV, QCD undergoes a chiral phase transition where 

5^ ■ the relevant degrees of freedom change from mesons and baryons to a phase of strongly coupled quarks and gluons, 

Oh the strongly coupled quark-gluon plasma (sQGP). This new phase of matter is being searched for in a number of 

■^r ' ultra-relativistic heavy ion facilities. There are a number of current observations in favor of the sQGP, ranging from 

1^ , hydrodynamical flow (soft probes) to jet quenching (hard probes). However, most of these observations are blurred 

ivj by the fact that the space-time evolution of the sQGP is short and its conversion to hadronic matter is involved. 

' Since the latter dominates the final stage of the evolution, it is producing competing signals that interfere with those 

y—i • from the sQGP. In this respect, dilepton and photon emissions are interesting probes of the collision region as neither 

J> I interact strongly with the medium produced in these collisions, thus they probe the early stages of the collision. This 

^^ is in contrast to hadronic observables which thermalize along with the collision region thus providing information only 

t^ on the late (or freeze-out) stage of the collision. 

Making quantitative predictions of the production rates of dileptons and photons is difficult for a number or reasons. 
Since the temperature produced in typical heavy-ion collisions is in the range of 200-300 MeV which is about the QCD 
^^ [ scale factor. A, the differential cross sections can not be computed in a weak-coupling expansion. Another uncertainty 
(^ ' is detailed knowledge of the evolution of both hadronic matter and quark gluon phase produced in heavy-ion collisions. 
^H ', In addition there is also a background of dileptons from other processes not occurring in the quark-gluon plasma such 
^_> . as hadronic decays. 

,_^ ' In the past there have been a number of experiments probing photons and dileptons created in hadronic collisions. 

O , One of the most recent experiments was the CERES (NA45) taking place as the CERN SPS collider which looked 

for dielectrons. It was found that the dielectron production exceeded the theoretical expectations for conventional 

processes in both hadronic and QGP matter Q, especially in the mass region 0.3 < M(GeV) < 0.6 p. A number 

^ ' of theoretical analyses were put forward to explain this excess based on effective Lagrangians with medium modifica- 

'^jf I tion [3, Q and dropping vector meson masses [3 ■ Model independent emission rates constrained by the strictures of 

''T^ . broken chiral symmetry and data were unable to account for the excess rate reported by NA45 [llH,!^. However, the 

C^ ■ large statistical and systematic errors reported by NA45 in exactly the excess region, did not allow for a definitive 

conclusion as to the theoretical nature of the emissivities. 

In this letter we revisit these issues in light of the recently reported dimuon data from the NA60 collaboration using 
In-In collisions at 158 Gev/Nuc |j|. These data have far better statistics, which gives much better constraints on any 
medium modification to the vector mesons [llL [l3| . We use the model independent analysis in Q, Q, |^ to analyze 
these data, whereby the emissivities are constrained by broken chiral symmetry in a dilute hadronic medium, and 
by non-perturbative QCD in the sQGP. The collision expansion and composition are extracted from an underlying 
hydrodynamical evolution set to reproduce the CERN SPS conditions. 

II. DILEPTON EMISSION RATES FROM A FIREBALL 

The rate of dilepton emission per unit four volume for particles in thermal equilibrium at a temperature T is related 
to the thermal expectation value of the electromagnetic current-current correlation function [l4L [13 . For massless 
leptons with momenta pi and p2, the rate per unit invariant momentum q = pi + P2 is given by: 

1^ - 3^Td?vr-w'(,) (1) 



where a — e^/47r, T is the temperature and 

W^(g) =1 f d'^x e*«-=^Tr[e^("^^«'^"^)/^r*J^(x)J,,(0)] (2) 

where eJ^ is the hadronic part of the electromagnetic current, H is the hadronic Hamiltonian, iib is the baryon 
chemical potential, N is the baryon number operator, and Q is the Gibbs energy. The trace is over a complete set of 
hadron states. 

In order to take into account leptons with mass rrii the right-hand side of Eq. ^ is multiplied by 

2171? , , 4to? , 1 /o 
(l + ^)(l-^)i/2 (3) 

To compare the theoretical dilepton production rates with those observed in heavy ion collisions, the rates must 
be integrated over the space-time history of the collision region and then finally integrated over the dilepton pair's 
transverse momentum and rapidity in order to compare with the yields measured by the NA60 collaboration. The 
final expression for the rates is given as: 

rIN f f f't" r°° /■''max /"Stt /-/4 d 

— = 2^M dy dq^ ■ q^ x Acc{M, q^,y) / rdr / drj rdr / de-^^-^(M, \q],T, ^b, x) (4) 

where M ~ ^ q^ is the dilepton invariant mass, y is the dilepton pair rapidity, r\ is the spatial rapidity, q\^ is 
the dilepton pair transverse momentum (with d defined as the angle between q\_ and the fluid element's velocity), 
X is the hadron fraction, r is radial coordinate (with r„iax set by the freeze-out temperature), and Acc{M^q\^^y) is 
the experimental acceptance taking into account that the CERES detector covers a limited rapidity in the interval 
y = 2.9 — 4.5 in the lab frame. 

The integration over yy, r, 9 and r was done over the full hydrodynamic simulation of the fireball as described below. 
\(f\ can be found by considering the two invariants q^q^ = A/^ and u^q^ constructed from the dilepton momentum 
and fluid 4-velocity which can be expressed as: 

q^ = (A/i cosh(y), q^,M^ sinh(y)) 

■u^ = (7_LCOsh(77),7^i;^,7xsinh(?7)) (5) 

giving 

l^l = [-Af 2 + (7^Afx cosh(?7) - u^q^ cos{e)f] ^'^ (6) 

where Ui_ = 7_lw_l, 7_l = , ^ , and w_l is the transverse fluid velocity which is taken from the hydrodynamic 

simulation. 

The acceptance function has a complicated dependence on Af , q±^ and y, but since our rates are y-independent we 
have used an acceptance with M and q±^ dependence built to speciflcations provided by the NA60 collaboration [l^ . 
Without detailed hadronic data available (such as transverse mass spectra and HBT analysis) a careful consideration 
of hadronic input, such as freeze-out temperature, cannot be made. Therefore there is a large uncertainty in the 
overall normalization of the yields, which depends strongly on Tfo.- In addition, the assumption of boost invariance 
can also affect the normalization as the acceptance probes very specific rapidities. The approach taken here is to 
normalize our results to the excess data in the peripheral centrality windows which fixes the normalization in the 
central bins. 

III. SPECTRUM ABOVE Tc 

At temperatures T>Tc lattice calculations have predicted that the relevant degrees of freedom consists of (strongly) 
interacting quarks and gluons. In order to compute the dilepton production rates as one would expect from a 
conventional phase of quark-gluon plasma we use the Born qq annihilation term [l9l |2(| which for massless quarks is 

q—u,d,s ' 
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FIG. 1: Integrated dimuon rates from the plasma phase for T=150 MeV. The thick sohd hne shows the perturbative qq 
annihilation rates while the this solid lines show the results for non- vanishing A|, B^ and E^ condensates and for only a non- 
vanishing A| condensate. For comparison the integrated hadronic rate at T=150 MeV and fj^B = 225 MeV is also shown, which 
will be discussed in the next section. 

where Nc is the number of colors, e^ the charge of the quarks, and n± = l/(e(*^l'?l)/^^ + !)• It should be mentioned 

that [7| reduces to the well-known vacuum result {ImW^ = ~'^^^ '^ ^ ei) at T=0 and that the finite temperature rate 
is always smaller then the T=0 rate due to Pauli blocking. 

It has also been seen in lattice simulations that near the critical temperature Tc there are still substantial chro- 
moelectric and chromomagnetic condensates present leading to additional non-perturbative effects. It was shown in 
|21| that the enhancement to the dilepton rates in a plasma with non-vanishing chromoelectric and chromomagnetic 
condensates can be given by 



ImW^^ = 



i(«c E «0[''<^^>4(t^V1(^b')1(^)(".(i-..)-.-(i-.-)) 



q—u,d,s 



(8) 



where the values of the above condensates in |H1 can only be estimated using non-perturbative calculations such as 
lattice QCD. The net result is a substantial enhancement (as seen in Fig. ^ for the case of dimuons) of the dilepton 
production rates below an invariant mass of ~ 0.4GeV. For the remainder of the paper we refer to the perturbative 
plasma rates as those given by H and the non-perturbative plasma rates as those given the sum of equations [3 and |HI 
using S = {^E^)/{200MeV)^ = {^B'^)/{200MeV)^ = and /? = (^^D/T^ = 0.4 which is the lower of the 
non-perturbative curves in Fig.^ An explanation of the choice of /3 = 0.4 can be found in |2ll |. 

IV. SPECTRUM BELOW Tc 

Even though there are various approaches to calculating production rates, they differ in the way in which the 
current-current correlation function in Eq. ^ is approximated and evaluated. The approach taken here is to use a 
chiral reduction formalism in order to reduce the current-current correlation function in |21 into a number of vacuum 
correlation functions which can be constrained to experimental e+e~ annihilation, r-decay, two-photon fusion reaction, 
and pion radiative decay experimental data. 

For temperatures, T < ttJtt and for nucleon densities, pN < Spo the trace in Eq. EJ can be expanded in pion and 
nucleon states. Keeping terms up to first order in meson and nucleon density gives j^] 



ImW^(g) = -3q^lmnviq^) + -^ f daWf{q,k) + f dNWf^ 



(q^p) 



(9) 



with phase space factors of 



dN = 



(Pp 1 



(27r)3 2Ep e(^p-'^^)/^ + 1 



(10) 



and 



da — 



d^k 1 



(27r)3 2uol e<IT^ - 1 



(11) 



with nucleon and meson energies of Ep ~ ^/ri'v^ + p^ and aj£ = \/'m\ + ^^ respectively. 

The first term in|51is the transverse part of the isovector correlator (0|T*VV|0) which can be determined experi- 
mentally from electroproduction data and gives a result analogous to the resonant gas model. At low and intermediate 
invariant mass the spectrum is dominated by the p(770MeT^) and /3'(1450MeT^). 

The term linear in meson density (the second term in Eq. ^^ can be related to experimentally measured quantities 
via the three flavor chiral reduction formulas 22] . It is shown in 0, Q that the dominant contribution comes solely 
from the part involving two-point correlators which gives: 



Wf (g, k) = i|g2lmn(,(g2) + ijg2j^n(n(.(q2) + Jn^(g2)) 



6 

T2 



6 



^(k + g)2lmn^((fc + qf) - ^(fc + qf\imYl\{{k + qf) + lmn^((fc + qf)\ + (9 ^ -q) 
4((fc • q? - mW)lu,Il\r(q^) X ^^^{k + q) + {q ^ -q) 



J 7T 



f 



((fc • g)2 - mlq^^ 



K 



)im(n(,(,^) + ^n^) 



x$RAj^(fc + g) + (g^-g) 



(12) 



Where A^ is the retarded meson propagator given by \/{q^ 



ie) and 11^ is the transverse part of the iso-axial 



correlator (0|T*jyiijyi|0). The spectral functions appearing in Eq. I|12l) can be related to both e+e" annihilation as 
well as r-decay data as was compiled in jl3i |. As already shown in JJ the spectral function can be directly related 
the form factor, Fy, via the KSFR relation where Fy is parameterized in the common Breit-Wigner form where the 
resonance parameters and decay constants are taken from empirical data. Included in the data are contributions to 
the spectral function from the p, uj, 0, oi, Ki and some of their radial excitations (see Table I in [3J). 

It can be seen in Fig.[21that the term linear in meson density decreases the rates from the resonance gas contribution 
for the mass region above the two pion threshold. However below the two pion threshold the only contribution to the 
rates come from the Ha terms in Eq. E| This is because the axial spectral density is integrated over all momentum 
in the thermal averaging (Eq. ^, which weakens the (fc + q)^ factor in Eq. 1121 allowing the 1/g^ term in Eq. QJto 
dominate at low q^. 

The final term in Eq. (Q which is proportional to the nucleon density is the spin-averaged forward Compton 
scattering amplitude of virtual photons off a nucleon. Experimentally, data is only available for values of g^ < 0, so 
while the photon rate which requires q^ — can be determined by use of the optical theorem the contribution to 
the dilepton rates must be determined by chiral constraints. Broken chiral symmetry dictates uniquely the form of 
the strong interaction Lagrangian (at tree level) for spin-i particles. Perturbative unitarity follows from an on-shell 
loop-expansion in -t' that enforces current conservation and crossing symmetry. To one- loop, the ttN contribution 
is parameter free. The large contribution of the A to the Compton amplitude near threshold is readily taken into 
account by adding it as a unitarized tree term to the one- loop result [3, Il0j | . The enhancement in the dimuon rates 
due to a non- vanishing baryon density can be seen in Fig. |21 where the solid curve shows the total dimuon spectra 
with an enhancement as large as a factor of two in the invariant mass region of 2m^ < Af (GeV) < 0.6. 



FIREBALL 



As mentioned earlier, in order to compare the theoretical dilepton production rates with those seen in heavy-ion 
collisions it is necessary to integrate these rates over the space time evolution of the collision region. We consider a 
region localized in space-time consisting of thermal hadronic matter acting as a source of particles. Equilibrium of 
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FIG. 2: The total integrated dimuon rates from a hadronic gas at T=150 MeV. The curve labeled Res. Gas shows the analogue 
of the resonant gas contribution (the first term in Eq. |UJ. The points labeled Tly and IIa give the contributions from all of 
the respective spectral functions in equations IHl and 1121 The thin line labeled meson is the total rate given by a hadronic gas 
without nucleons. The thick solid line gives the total dimuon rate when nucleons (shown here for hb = 225 MeV) are taken 
into account. 

the collision region is strictly a local property with different temperatures and baryon densities possible in different 
space-time domains. 

A computational hydrodynamic code was already developed by one of us and it has been modified to the conditions 
of the SPS collider for fndium on Indium collisions. In this paper we only briefly outline the physics behind this code 
and show the results of the In-In collisions which has not been modeled before. For technical details regarding the 
hydrodynamic calculations the reader is referred to the prior works by one of us |2J| . 

A. Hydrodynamics 

The hydrodynamic equations for a relativistic fluid consist of the local conservation of energy and momentum, 
which can be written in compact form as 9„T'"' — 0, as well as local charge conservation 9^ Jf — where T^" — 
(e + p)U^U'^ — pg^'^ is the energy-momentum tensor with e the energy density, p the pressure, f/^ — 7(1, v) is the 
proper velocity of the fluid, and Jf is any conserved current [e.g. isospin, strangeness and baryon number in the case 
of strong interactions). 

The same space-time evolution scenario as first proposed by Bjorken ^S] is assumed where the equation of motion 
can be described by the Bjorken proper time t = ^JW—z^ and the spatial rapidity y = -^Awj^. One of the main 
results, following from the assumption of a central-plateau structure in the rapidity distribution is that of boost 
invariance, stating that the initial conditions and thus the subsequent evolution of the system are invariant with 
respect to a Lorentz boost. Thus a solution at any value of y can be found by boosting the solution at y = to a new 
frame moving with velocity v = tanh(?/) in the negative z-direction. 

With the assumption of boost invariance the equations of motion are a function of the transverse coordinates and 
the proper time r only. After integrating over the transverse plane of the collision region one finds that {dStot/d-q), 
{dns/dri), and the net transverse momentum per unit rapidity are all conserved. 



B. Equation of State 



In order to solve the equations of motion as given by the vanishing of the divergence of the energy-momentum stress 
tensor one must have an Equation of State (EoS) relating the local values of the pressure, energy density, and baryon 
density (ns). The approach taken here is to consider an EoS with a variable latent heat in the e/n^ plane where the 



Parameter 


Value 


'1 


0.05c 


^QGP 


0.33c 


Tc 


170 MeV 


T/.o. 


130 MeV 


TO 


1.0 fm/c 


ub/s 


0.0238 


Cs 


8.06 


(-■ng 


0.191 



TABLE I: Parameters used in the hydrodynamic simulation of In-In collisions, 
following two derivatives hold along a path of constant ns/s: 

(t) = -^ (») 

If the speed of sound is defined everywhere along with the entropy of one arc in the e, ub plane the above derivatives 
can be integrated in order to determine s(e, ub)- From the entropy all other thermodynamic variables, such as T and 
fiB, can be found as needed. 

We consider a fireball that consists of three phases, a hadronic phase, a QGP phase, and a mixed phase. The 
hadronic phase is taken to be made of ideal gas mixtures of the lowest SU(3) multiplets of mesons and baryons. All 
intensive thermodynamic quantities including p, e, s, and ub can be found as a sum of that quantity's contribution 
from each specie in the gas consisting of a simple Bose or Fermi distribution. The hadronic phase is taken up to a 
temperature of Tc < 170 MeV and has a squared speed of sound of approximately l/5c^. For temperatures above 
Tc only the squared speed of sound, c^ is specified. For the mixed phase it is taken almost at zero (c^ = 0.05c). For 
the QGP phase the degrees of freedom are taken to be massless and the speed of sound is accordingly c^ = l/3c. The 
extension of this analysis to the sQGP is beyond the scope of this work. 

C. Initial Conditions 

The initial conditions of the fireball consist of setting the entropy and baryon density proportional to the number 
of participating nucleons in the transverse plane at some initial proper time tq = 1 fm/c. Since both the entropy and 
baryon number per unit rapidity are conserved the final yields of pions and nucleons are proportional to the number 
of participants. The number of participants were calculated by use of a Glauber model and the initial entropy and 
baryon densities were fixed by two constants Cs and Cng , which respectively are the entropy and net baryon number 
produced per unit spatial rapidity per participant. These constants were fixed to the conditions at the CERN SPS 
collider in order to fit the total yield of charged particles and the net yield of protons. Tabic Q] summarizes the input 
parameters used in the hydrodynamic calculations. In order to address the centrality of the collision the impact 
parameter was chosen in order to reproduce the number of participants as reported in |25j . 

D. Results for In-In Collisions at CERN SPS 

The hydrodynamic result for In-In Semi-Central collisions is shown in Fig. O The two thick lines labeled eg and 
en represent contours of constant energy density showing the transition from the plasma phase to the mixed plasma 
and hadronic phase and the transition from the mixed to the purely hadronic phase respectively. It can be seen that 
the QGP phase takes up a much smaller space-time volume then the hadronic phase, however the rates still appear 
in the spectrum as the high temperatures in this region enhance the rates by an order of magnitude. The effect of 
nucleons depends on the baryon chemical potential in the fireball. This is plotted as a function of temperature for 
the pure hadronic phase in Fig. ^ 
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FIG. 3: The hydrodynamic solution for semi-central In-In collisions at the SPS collider. The thin lines show contours of constant 
transverse fluid rapidity {v± — tanh{y±)) with values of 0.1, .02, ..,0.5. The dashed lines show contours of constant temperature 
with values of (working radially outward) T=200 MeV, T=150 MeV and T=135 Mev. The eg = 1.70 GeV/fm^ and en = 0.50 
GeV/fm'^ contours represent the phase changes from QGP to mixed and from mixed to hadronic matter respectively. 




FIG. 4: Dependence of the baryon chemical potential, /is, on the temperature for the hadronic phase of the fireball. 



VI. RESULTS AND DISCUSSION 



Our final diinuon yields for the four centrality windows is shown in Fig. I VII where it is compared to the excess data 
measured by the NA60 collaboration. In all four figures we show the total dimuon yield, which includes contributions 
from the hadronic phase, either the perturbative or non-perturbative plasma phase, as well as the DD contribution 
as provided by the NA60 collaboration. For all cases we also show separately the perturbative and non-perturbative 
plasma contributions to the overall yield. It can be seen to be almost negligible in the peripheral data. For the central 
data, where there is a larger plasma contribution we also show curves showing separately the hadronic contribution. 

Even though it can be seen that the theoretical rates are able to describe most of the features of the spectrum, a 
number of things should be noted before a direct comparison is made. The rates below M=0.4 GeV should not be 
taken literally since they are obtained by saturating the total measured yield in that region by ij Dalitz decays only, 



8 



thereby lowering the excess close to the dimuon threshold. Actually, by reducing the rj yield by 10% the data has 
much better agreement with the theory for M< 0.4 GeV. The charm decay data was analyzed and provided by NA60, 
and since the contribution from charm decays is not subtracted from the excess data it must added to our rates for 
comparison with experiment. The excess spectra which is shown in the figure is created from subtracting the cocktail 
(omitting the p(770)) from the total observed data. This would erase any uj or (f) peak at the vacuum positions. Since 
our hadronic rates don't modify either the position or width of the uj or (j) it can be very difficult to distinguish any 
residual ui or 0's from the cocktails'. 

It can be seen right away that the dimuon yields are reproduced in the peripheral centrality windows. This is 
expected as the matter is dilute and any medium modification to the spectral densities will be accounted for in the 
virial expansion (Eq. O . In the central bins it can be seen that the shape of the spectrum changes as one goes to 
more central collisions. Even though the general shape of the spectrum is reproduced by our rates, our rates slightly 
over-predict the yield at the p peak by about « 50% for semi-central and by « 60% for the central data. Even though 
our rates agree fairly well with the remaining data away from the p peak, there is still room for enhancement in the 
low mass region, 0.4 < M (GeV) < 0.6. 

We should finally mention what happens when the non-pcrturbative QGP rate is used instead of the perturbative 
result. Similar to the perturbative QGP results in Fig. IVII the non-perturbative plasma rate is about a factor or 
two larger in the low mass region. Even though this does help to explain some of the excess in the low mass region, 
especially in the more central data, the space time volume of the plasma phase is too small to have a large effect. 
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FIG. 5: NA60's excess dimuon data compared to our thermal yields which include contributions from either the perturbative 
or non-perturbative QGP phase, the hadronic phase and the DD contribution. Shown for all four centrality windows. 



VII. CONCLUSIONS 

Using a parameterization of the results given by a hydrodynamic model of the collision region at the CERN 
SPS collider, the NA60 dimuon spectrum was reproduced using a pure thermal model assuming that there exists a 
sQGP phase above Tc with an interacting hadronic phase persisting until freeze-out. The dimuon spectrum from 
the sQGP phase originates primarily from qq annihilation with non-perturbative effects due to non-vanishing gluon 
condensates. After hadronization it is assumed that there remains a dilute hadronic gas in which the dimuon rates 
can be constrained entirely from broken chiral symmetry arguments and experimental data. The combination of these 
two rate equations, after being folded over the space-time evolution of a fireball, are able to explain most of the excess 
dimuon data, especially in the more peripheral collisions where our assumptions about diluteness hold. In the more 
central data, where the assumption of diluteness may breakdown, it is necessary to investigate how higher order terms 
in the virial expansion modify the spectrum. 
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